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Abstract 



We report on our observations of the gamma-ray and optical photopolari- 
metric behavior of the radio-loud, narrow line type-1 Seyfert galaxy PMN 
J0948+0022 over a sixteen-month period. As this object has recently been sug- 
gested to represent a prototype of an emerging class of blazar-like objects, the 
observed properties are compared to those of blazars. While our observations 
do suggest a correlation between the gamma-ray and optical behavior, they do 
not offer conclusive proof that the emissive behavior of PMN J0948+0022 is due 
to a relativistic jet oriented close to our line of sight. We thus advocate for, 
and describe the requirements of, continued monitoring of this source in order to 
better address this issue. 

Subject headings: galaxies: active - galaxies: individual: J0948+0022 - galaxies: 
photometry - galaxies: polarimetry - galaxies: Seyfert 
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Introduction 



Recently, a new sub-class of active galactic nuclei have been observed with properties 
that would have previously been divided among Seyfert galaxies and Broad Line Radio 
Galaxies (BLRG), such as blazars. These objects, the Radio-Loud Narrow Line Type-1 
Seyferts (RL NLS1), posse ss the standard identifying properties of NLS1 identified by 



Osterbrock k Poggel (119851 ): FWHM(Hp) < 2000 km/s, strong optical emission of Fell, 
and weak emission from forbidden lines (i.e. [QUI] /Hp < 3). However, the property of 



radio-loudness (R > 10, where R = fi.4 ghz/^ooa) (IKellermann et al 



1989) is markedly 



rare in galaxies of this type, occurring in only < 7% of such systems (IKomossa et al. 
High brightness temperature (> 10 13 K) as well as loudness in the radio ( jZhou et al. 



2006). 



20031 ) 



as well as shorter wavelengths up to and including y-rays, and strong/rapid variability 
across short (minutes) to long (years) time scales are, however, properties of blazars. It is 
the combination of these various observational properties that lead many to now believe 
that RL NLSyl and blazars may both play hosts to relativistic jets. 

It is now widely accepted that all varieties of Active Galactic Nuclei (AGN) are 
manifestations of the same basic phenomenon - accretion of matter onto a Supermassive 
Black Hole (SMBH) at the center of a galaxy. The different classes of AGN that we then 
observe result, to a large degree, from these objects being oriented differently with respect 
to our line of sight. Blazars, a class of AGN characterized by strong and variable emission 
across all wavelengths an d strong and highly variable polarization in the radio and optical 



( IBlandford &: Rees 



19781 ). are believed to result from our line of sight near the axis of 



a relativistic jet of particles being emitted from the central engine of the source. Until 
recently, blazers were almost exclusively o bserved to be hosted in elliptical galaxies, with 



very few exceptions (IMcHardy et al 



19941 ). As the majority of NLSyl hosts are spiral 



galaxies, finding evidence of blazar-like behavior in such systems would help to fill in a 
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curiously barren demographic of the blazar population. 



PMN J0948+0022 is one such object that displays these varied properties. It possesses 



the e xpected properties of a Narrow-Line Seyfert-1 galaxy as described above ( iZhou et al. 



20031 ). as well as those of blazars, such as 



strong and variable emission in t 



gamma-ray energies over long timescales (lAbdo et a. 



microvariability in the optical (IMaune et al 



20131 ) 



20091 



Foschini et al- 



ike] iri et al 



re radio thru 



20121), and 



(1201 ll ) also observed 



J0948+0022 to exhibit a very high degree of linear polarization (18.8%) in the optical 
(V-band). However, their measurement was obtained while the object underwent a .74 
magnitude change in brightness in a single observing sequence, thus there is no way 
currently to corroborate this impressive result. 

To date, a comprehensive, long-term study of the optical polarimetric characteristics 
of RL NLSyl has yet to be done: that is one focus of this work. By combining 
optical polarimetric observations with data from other regimes (namely, high-to-medium 
cadence optical photometry and gamma-ray photometry from FERMI/LAT), we will 
be able to obtain a more comprehensive picture of the nature of this emergent class of 
relativistically-jetted AGN. 



2. Observations and Data Reductions 

2.1. Optical Photopolarimetry Data 

All optical photopolarimetric data used in this study were obtained with the 72" 
Perkins telescope at Lowell Observatory in Flagstaff, Arizona, using the PRISM instrument 
which includes a polarimeter with a rotating half-wave plate. Data were obtained during 
seven observing runs between February, 2011 and May, 2012. The specific dates of each 
observation are given in Table 1. The observations consisted of a series of 2-4 measurements 
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for the Q and U Stokes parameters per object. Each series consisted of four images, each 
taken at different instrumental position angles - 0°, 45°, 90°, and 135° - of the waveplate. 

Corrections to polarimetric values were obtained from two sources: in-field comparison 
stars and seperately- observed polarimetric standards, both polarized and unpolarized 



( iSchmidt et al. 



19921 ). As the camera has a wide field of view (approx. 14' x 14'), we are 
able to use field stars for interstellar polarization corrections by subtracting the average 
percent polarization of the brightest field stars. Polarized and unpolarized standard stars 
are used to calibrate corrections for polarization Position Angle (P. A.) and instrumental 



polarization (typically less than 1%) (IJorstad et al. 



20101 ). respectively. 



The data were reduced and analyzed using in-house scripts, which utilize standard 
packages in the PyRAF 2.0 suite of reduction toolj*]- Bias frames were taken at the 
beginning of every night and combined into a master bias that was subtracted from each 
image. Flat frames were taken at least once per run, using a featureless screen inside the 
dome. Each position of the waveplate required its own set of flats, which would later be 
combined into one master flat per position angle for application to the appropriate science 
image(s). Cosmic ray cleaning was performed on all science images, with the threshold 
and fiuxratio parameters set to 35 and 5, respectively. Aperture photometry was then 
performed on the calibrated science frames on an object-by-object basis. An aperture radius 
of 7 arcsec was used on all images both to maximize the signal-to-noise an d to maintain 
consi stency with the optical photometry being performed on this target by 



Mauneetal 



( 120131 ) . Use of the in-field comparison stars compiled by Maune et al. allowed for the ability 
to obtain simultaneous measures of the absolute R-band magnitude, as well as the percent 
polarization (P) and electric vector position angle (EVPA). 



x PyRAF is a product of the Space Telescope Science Institute, which is operated by 
AURA for NASA. 
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For ease-of-comparison with gamma-ray data, optical data was converted from 
magnitudes to units of flux (MeV) using the following equation: 

F = 2941 * io-°- 4 * Ma s 

Where F is the flux in mJy and Mag is the R-band magnitude. Fluxes in mJy were 
then converted to MeV for direct comparisons to gamma-ray data. 



Maune et al 



(120131 ) was also included in 



R-band photometry which we presented in 
this study. As the optical data were also obtained at a significantly higher (as well as more 
variable) sampling cadence than the gamma-ray data, all optical data taken during a single 
observing run were averaged into a single data point. The total optical lightcurve presented 
herein thus consists of 28 data points, covering a time period from February 10, 2011 to 
May 19, 2012. 



2.2. FERMI-L AT Data 

Gamma-ray data were obtained through the FERMI-LAT public data server. The 
Large Area Telescope (LAT), on board the Fermi Gamma-ray Space Telescope, is a 
pair-conversio n detector sensitive to gamma-rays in the 20 MeV to several hundred GeV 



energy range ( lAtwood et al 



20091 ). The instrument has worked almost continuously in 
all-sky-survey mode since its launch in 2009, which allows coverage of the entire gamma-ray 
sky approximately every 3 hours. The data were reduced and analyzed using Science Tools 
v9r27pl and instrument response functions P7SOURCE_V6. We utilized the likelihood 
analysis procedure as described at the FSSC website. Photon fluxes in Table 2 are calculated 
using data from MJD 55562 to MJD 56179 (Jan. 01, 2011 to Sep. 08, 2012). 

Our data were downloaded from the FERMI website on September 8, 2012 and cover 
a region of interest (ROI) on the sky 15° in radius, centered on the location of PMN 
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J0948+0022 (2FGL0948.8+0020 from the Fermi 2- Year Point Source Catalog), and in an 
energy range of 100 MeV to 300 GeV. Our gamma-ray light curve consists of 20 equally-sized 
bins, each of which is 29.530589 days in length so as to be centered on the new moon, 
which was the approximate time that our observing runs at Lowell Observatory took place. 
The first bin began on February 03, 2011, while the last bin ended on September 8, 2012. 
Only data corresponding to the SOURCE class (evclass=2) was utilized, with a 52° cut-off 
rock-angle of the spacecraft, while an additional cut utilizing an angle of 100° from the 
zenith was imposed so as to minimize the contamination due to gamma-rays coming from 



Earth's upper atmosphere. As PMN J0948+0022 is within 20° o 



the ecliptic, and the Sun 



201 ll ) . a final cut was used 



is a source of gamma-rays comparable to our target (lAbdo et al.l 
to exclude exposures that occurred when the Sun was within the 15° ROI. Photon fluxes 
and spectral fits were derived using an unbinned maximum likelihood analysis which was 
accomplished using the Science Tool GTLIKE. 

In order to accurately measure the flux and spectral parameters of the source, 
gamma-rays emitted from the background needed to be accounted for. To this end, two 
models were used: and isotropic background model accounting for extragalactic diffuse 
emission and residual charged particle background, and a Galactic diffuse emission model 
to account for diffuse sources from within our own galaxy. The isotropic model we used was 
the one contained in the file iso_p7v6source.tx1o, while the Galactic component was given 
by the file gal_2year7v6_v0.fits. The normalizations of both components were left to vary 
freely during likelihood analysis. 

In order to determine the significance of the gamma-ray signal from PMN J0948+0022, 
we used the Test Statistic (TS). The Test Statistic is defined as T S _ 2Alog(likelihoo d). 



where likelihood refers to the likelihood ratio test as described in iMattox et al. 



1996 



2 http: / / fermi.gsfc.nasa.gov / ssc / data / access /lat /Background /Mo dels, html 
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Determining the likelihood of a given photon flux being produced by a source with a 
given spectral model was accomplished using the GTLIKE Science Tool. Our source 
model consisted of all the known gamma-ray point sources located within a 15° radius 
of 2FGLJ0948. 8+0020. Initial values for all spectral parameters for these sources were 
taken from the LAT 2-year Point Source Catalog. Along with PMN J0948+0022 and the 
aforementioned background models, several point sources were allowed to vary (i.e. photon 
indices, normalization factors, and spectral slope indices were left as free parameters) 
during the likelihood analysis, so as to account for the inherent variability of many gamma 
ray sources. The type of spectral model u sed for a given sou rce was the same model used 



for that source in the LAT 2-year catalog ( jNolan et al. 



20121 ). 



For PMN J0948+0022 specifically, we used a LogParabola model to describe the 
gamma-ray spectrum of the source in this study. This model takes the form: 

N{E) = AT (^)"+Pln(E/E i ,) 

where No is the normalization index, a is the photon index of the pivot energy Eb, and 
P is the curvature index. All of these parameters were left un-fixed with the exception of 
Eb: this was fixed at a value of 271.597 MeV, which is the value reported for this object in 
the LAT 2-year catalog. All gamma-ray fluxes presented in this paper are integrated over 
the entire energy range cited above (100 MeV to 300 GeV), with corresponding units of [ph 
cm~ 2 sec _1 ]. 

Sources outside our 15° Radius of Interest (Rol) but within 25° of the target were also 
included in the source model, as the point spread functions of these objects could result 
in extra photons seeping into the Rol of our target. All parameters for these sources were 
fixed to their 2FGL catalog values during the analysis. 
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3. Results 

Our polarimetric results are detailed in Table [I] and displayed for comparison in Fig. |2j 
Note that all photometric data points in the aforementioned table and figure were derived 
from polarimetric measurements. PMN J0948+0022 displays a moderate yet significant 
amount of variability in terms of percent polarization and EVPA, with a maximum P of 
8.17 ± .74% and an EVPA which varied widely. No significant correlation was observed 
between the optical or gamma-ray fluxes and the polarimetric quantities, though there does 
appear to be some correlation between changes in the percent polarization and EVPA. 
However, it should be noted that the 180° uncertainty inherent in the measurement of the 
position angle may produce the appearance of trends where none exist. 

The optical photometric data obtained for this study are also presented in Table [IJ 



Maune et al. 



(12013( 1 . 



while additional optical data used in this study was presented in 
The two datasets were merged to make a master optical light curve, which contained 953 
R-band data points. This dataset consisted of data from several observing runs, five of 
which accounted for over 90% of the data. Therefore, a program was written which created 
bins from the data in which no two adjacent data points were separated by more than one 
day's time. An average value was then calculated for each bin to provide a single data 
point (28 optical photometric data points in total), allowing us to better- match the optical 
photometric data to the cadence of the FERMI observations. 

The long-term behavior of PMN J0948+0022 in gamma-rays is also shown in Fig. [TJ 
Table [2] tabulates the gamma-ray data, showing the date, TS value, and photon flux of 
each data point. It is common practice to use the TS value as a proxy for the confidence in 
measurement, with TS > 25 (corresponding roughly to a > 5) being the threshold below 
which an upper limit is calculated. We applied the method for extracting upper limits, 
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as described on the FERMI websitd_|, to the two data p oints with t 



le lowest TS value. 



20091 ). and only two 



As PMN J0948+0022 is a well-established y-ray source flAbdo et all 
of our data points are significantly below this threshold, we treat all of our remaining 
gamma-ray data points as direct detections. It should also be noted that the exclusion of 
exposures which contained both the Sun and PMN J0948+0022 had very little effect on 
our results, with the exception of the bin centered on MJD = 5802.21 (August 28, 2012), 
whose measured flux was virtually unchanged while the associated TS value was reduced by 
roughly a factor of five. 

An essentially simultaneous rise was observed for this object in the optical R-band 
and in the gamma-rays, followed by a return to a faint state as shown in Fig. [TJ PMN 
J0948+0022 was observed to undergo a significant increase in gamma-ray brightness in the 
spring of 2011, which was mirrored by a steady increase of the average R-band luminosity 
over the next few months. Unfortunately, the gamma-ray flare which was observed, 
continued into the summer of 2011 when PMN J0948+0022 could not be viewed in the 
optical due to the close proximity of the Sun. Nevertheless, when optical monitoring 
resumed in November 2011, the object was observed to be back at a low state in both the 
optical and gamma-rays. 



4. Conclusions 



Mu lti-epoch VLBA ob servations of a number of gamma-ray detected quasars and 



blazars ( Uorstad et al. 



120011 ) suggest that the gamma-ray emission observed for these objects 
originates near the radio core, perhaps corresponding to a standing shock in the jet itself, 
and not specifically originating from a location near the central SMBH. It is thought that 



? http: / /fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/python_tutorial. html 
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the oft-observed flares are produced as a result of turbulence/instabilities, which are present 
in the relativistic jet and arise as shocks in this outflow. With the passage of a shock down 
the jet, the magnetic field will be compressed within the shock region causing the field to 
become more highly ordered. Observationally, one would then expect that the fractional 
linear polarization would increase and the direction of the electric vector position angle 
would change rapidly. 

The present observations allow us to evaluate if this is the case for the gamma- 
ray/optical flux variations observed for PMN J0948+0022. In Fig. [T]we have plotted the 
gamma-ray and optical light curves for 2011-2012. We see, visually, that the optical and 
gamma-ray variations seem loosely related. In order to evaluate this, in Fig. [3] we have 
plotted the optical flux versus the gamma-ray flux. As is easily seen in this figure, there 
is only a weak relationship between the gamma-ray and optical flux during the time of 
these observations. The absence of a strong correlation suggests that there is substantial 
turbulence present in the jet and the magnetic field is not highly ordered, suggesting that no 
strong shock is present during the time of these observations. This is supported by the data 
in Table [U which shows that at the time of the brightest optical observation (May 24, 2011) 
there was a very low polarization (P=1.35%). The highest polarization was observed in 
March 21, 2012 (P=8.17%), with significant changes in both P and EVPA occurring within 
a few days before and after this occurrence. However, no significant flare was observed in 
either optical or gamma-ray flux near this time. 

While the present observations do not allow one to definitively confront the model 
suggesting that the gamma-ray emission is produced at a shock, down-stream in the 
jet, some distance from the SMBH, it does suggests what is required: It will require 
quasi-simultaneous gamma-ray and optical/photopolarimetric observations during a major 
gamma-ray/optical flare consisting of a change in flux significantly greater than a factor of 
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three, such as was observed in the present campaign. Under these conditions, one should be 
able to determine if there is a significant increase in the polarization (P) and the expected 
rapid change in the electric vector position angle (EVPA) accompanying such an outburst. 
Therefore we encourage continued photopolarmetric monitoring of this object in order 
to investigate the behavior of the polarization during the next major outburst of PMN 
J0948+0022. 



The authors would like to thank Svetlana Jorstad and Paul Smith for offering their 
their assistance and technical expertise. 
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PMN J0948 + 0022: 2011-2012 R-band St 7- ray Data 

1 1 1 1 1 '■ r 




5600 5700 5800 5900 

Time (MJD) 
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6100 



6200 



Fig. 1. — Gamma-ray flux (integrated from lOOMeV to 300 GeV) over-plotted with averaged 
R-band flux (blue dots and green x's, respectively), downward-pointing arrows denote upper- 
limits on gamma-ray data points. 
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Date 


R-mag (err) 


P ferr") 


EVPA (err) 


PpVi 1D 9D1 1 
reu. iUj zuii 




9 09 (r\ ^rh 


^fi (^ q 99"i 


Mav 94 901 1 


18 91 1 (0 0191 


1 35 fl 99") 


93 49 (R fi1 1 


Jan. 25, 2012 


19.053 (0.027) 


1.66 (0.26) 


-50.74 (2.6) 


Feb. 25, 2012 


18.677 (0.030) 


1.70 (0.48) 


19.37 (13.64) 


Feb. 26, 2012 


18.824 (0.028) 


5.95 (0.50) 


-59.27 (1.98) 


Feb. 27, 2012 


18.613 (0.033) 


2.45 (0.16) 


19.93 (5.26) 


Mar. 21, 2012 


18.709 (0.026) 


8.17 (0.74) 


58.35 (13.52) 


Mar. 22, 2012 


19.479 (0.024) 


1.92 (0.39) 


-46.67 (0.26) 


Apr. 22, 2012 


19.065 (0.024) 


3.89 (0.61) 


28.94 (13.48) 


Apr. 23, 2012 


18.887 (0.027) 


3.13 (0.36) 


8.99 (1.71) 


May 12, 2012 


19.191 (0.027) 


0.90 (0.28) 


6.20 (62.69) 



Table 1: Photopolarimetric observations of J0948+0022 obtained between February, 2011 
and May, 2012. Columns are: (1) time of the observation in MJD (JD - 2.45e6), (2) optical 
R-band magnitude and (error), (3) percent polarization of target and (error), and (4) EVPA 
and (error). 
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Date 


TS 


Flux (err) 


Date 


TS 


Flux (err) 


5595.50 


21.44 


0.622 (0.260) 


5890.81 


29.48 


0.787 (0.230) 


5625.03 


11.85 


0.789* 


5920.34 


35.16 


1.160 (0.033) 


5654.56 


55.03 


1.374 (0.356) 


5949.87 


38.95 


1.275 (0.306) 


5684.09 


129.86 


2.403 (0.353) 


5979.40 


25.52 


0.697 (0.306) 


5713.62 


59.15 


1.846 (0.382) 


6008.93 


41.26 


1.046 (0.286) 


5743.15 


201.94 


2.471 (0.317) 


6038.46 


13.86 


0.693* 


5772.68 


198.17 


2.847 (0.387) 


6067.99 


60.03 


1.135 (0.125) 


5802.21 


21.34 


1.633 (0.500) 


6097.52 


100.94 


2.083 (0.351) 


5831.74 


189.07 


2.271 (0.314) 


6127.05 


82.85 


1.610 (0.329) 


5861.28 


213.34 


2.853 (0.058) 


6156.58 


23.10 


1.463 (0.514) 



Table 2: The results of the likelihood analysis for each bin. Columns are: (1) the midpoint 
of the observation (MJD), (2) the TS value of the measurement, and (3) the integrated flux 
from 100 MeV to 300 GeV for each bin and the (uncertainty) on that value. Values with a 
* are upper limits only. 
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PMN J094S+0022: 2011-2012 Optical, Pol, & 7 -ray Data 
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Fig. 2. — A comparison of R-band optical data obtained with polarimetry (top panel), 
the Percent Polarization (second panel), position of the electric vector in degrees (third 
panel), and the integrated gamma-ray flux (bottom panel) of PMN J0948+0022. The same 
horizontal axis is common to all four plots. 
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PMN J094S+0022: 2011-2012 R-band vs. 7 -ray Flux 
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Fig. 3. — A plot of concurrent R-band and linearly-interpolated gamma-ray fluxes. 



